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New Self-Similar Solutions of Polytropic Gas Dynamics 
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ABSTRACT 

We explore semi-complete self-similar solutions for the polytropic gas dynamics in- 
volving self-gravity under spherical symmetry, examine behaviours of the sonic crit- 
ical curve, and present new asymptotic collapse solutions that describe 'quasi-static' 
asymptotic behaviours at small radii and large times. These new 'quasi-static' solutions 
with divergent mass density approaching the core can have self-similar oscillations. 
Earlier known solutions are summarized. Various semi-complete self-similar solutions 
involving such novel asymptotic solutions are constructed, either with or without a 
shock. In contexts of stellar core collapse and supernova explosion, a hydrodynamic 
model of a rebound shock initiated around the stellar degenerate core of a massive pro- 
genitor star is presented. With this dynamic model framework, we attempt to relate 
progenitor stars and the corresponding remnant compact stars: neutron stars, black 
holes, and white dwarfs. 

Key words: black hole physics — hydrodynamics — shock waves — stars: neutron 
— supernovae — white dwarfs 



, 1 INTRODUCTION 



X 



The self-similar evolution of gas dynamics with spherical 
symmetry under self-gravity and thermal pressure has 
attracted considerable i nteres t in con te xts of various astro- 
physic al applications [Shu] |197] J ): | Goldreich fc Weberl 
<198dk iFillmore fc Goldreichl Y\98$i : |Lou| J2005D : 
ll^i^rierrLou^^Zh^n^ni20uf3) l. There are several rea- 
sons to pursue such model studies: numerical simulations 
show the tendency of self-similar behaviours in spherical 
flows (e.g. Bodenheimer & Sweigart 1968 and Foster & 
Chevalier 1993); this model framework is sufficiently simple, 
yet can still carry valuable information on the roles of the 
competing forces of gravity and thermal pressure. Larson 
(1969a, b) and Penston (1969a, b) were the first to analyse 
such similarity problems in astrophysical contexts. Since 
then, different asymptotic solutions and global numerical 
soluti ons have be en found and developed. In the isothermal 
case, IShul il977l) obtained the inner free-fall asymptotic 
solution with divergent speed and density profiles at small 
x, and constructed the expansion wave collapse solution 
(EWCS) to suggest the inside-out collapse scenario f or the 
formation of low-mass stars. Subsequently, iHunterl jl977l) 
discovered a sequence of discrete solutions, referred to as 



* E-mail: louyq@tsingliua.edu.cn and lou@oddjob.uchicago.edu; 
wwg03@mails.tsinghua.edu.cn 



the Hunter type solutions, which have finite density and 
speed profiles at small radii, by a matching procedure in 
the d ensity-speed phase diagram. IWhitworth fc Summers! 
(1985) noted the existence of asymptotic solutions of 
constant speed at large radii, as a generalization of the 
asymptotic solution at large x of Shu (1977), and con- 
structed solutions with weak discontinuities across the 
sonic critical line. The propertie s of such wea k disconti- 
nuities were f urther discussed by IHunterl l|l986i) . Recently, 
iLou fc Shenl J2004T) constructed isothermal solutions to 
describe self-similar evolution of envelope expansion with 
co re collap s e (EE CC) using the similar matching procedure 
of IHunterl Jl977l) . Bian & Lou (2005) explored various 
similarity isothermal shock solutions. 

For a polytropic gas, IChend lll97ct) introduced a gener- 
alized self-similar transformation in the sense of n = pjp 1 
being constant along each streamline but not globally, where 
1^7^5/3 and the initial density profile may be ad- 
justed, and mainly studied the polytropic counterparts of 
isothermal EWCSs. lYahill (1983) developed a model for n 
being a global constant, discussed the polytropic case for 
6/5 ^ 7 ^ 4/3, and mainly focussed on the pol ytropic coun- 
terpar ts of Hunter-type isothermal solutions. ISuto fc Silld 
( 1988) used a similarity transformation similar to that of 
Yahil dl983l) . but discussed a generalization of the equation 
of state for 1 < 7 < 4/3 involving yet another parameter n, 
and mainly constructed the counterparts of isothermal free- 
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fall solutio ns both crossing and no t cross ing the sonic crit- 
ical curve; iMcLauehlin fc Pudrit j il997l) treated the limit 
case of 7 — *• using the logotropic equation of state, and 
considered the free-fall s olutions and EWCSs as c o unter - 
parts of earlier analyses. iFatuzzo. Adams fc Mversl <l2004l) 
explored the possibility of an initial equation of state being 
different from a later dynamic equation of state, and exam- 
ined the dependence of the accreted mass versus the poly- 
tropic index; their self-similar solutions do not encounter 
the sonic critical curve. The basic polytropic model frame- 
work has a lso been extended in var i ous re levant physical as- 
pects, e.g.. lTere bev, Shu fc Cassenl dl984f) for a slowly rotat- 
ing gas cloud. iBoilv&^Lvnderi^B elll (^?9 5p for including var- 
ious f orms of radiative losse^nSemeh^ ^Mcjiez fc de Vegal 
(2001) for including viscosity. IWan^^^LmTil200^) ^r mod- 
elling a random magnetic field, as well as other relevant 
research works involving hydrody namic and MHD shocks 
jBian fc Loull2005l:lYu et al]l20od) . 

Among these earlier research works. ISuto fc Silkl (Il98ct) 
introduced a straightforward self-similar transformation and 
analyzed polytropic self-similar flows as a generalization of 
an isothermal gas (e.g. Shu 1977). They analyzed both cases 
of 7 > 1, n = 1 and of n — 2 — 7. The n = 1 case is 
somewhat artificial because this equation of state evolves 
with time and does not conserve the local specific entropy. 
We therefore focus on the n = 2 — 7 case in this paper 
to explore a polytropic model for the collapse of a mas- 
sive star. As shocks are ubiquitous and important in vari- 
ous astrophysical processes (e.g. Kennel & Coroniti 1984 for 
an MHD model of the Crab Nebula with shocks), and self- 
similar shocks have been investigated in various astrophys- 
ical contexts (e.g. Tsai & Hsu 1995; Shu et al. 2002; Shen 
& Lou 2004; Bian & Lou 2005; Yu, Lou, Bian & Wu 2006), 
we further construct self-similar shocks in a polytropic gas 
in this paper. We analyze the behaviour of the sonic crit- 
ical curve and present new asymptotic solution referred to 
as the 'quasi-static' asymptotic solution. Semi-complete so- 
lutions containing the 'quasi-static' asymptotic solution as 
the limit are constructed, either with or without a shock. We 
invoke this simple theoretical framework to model the col- 
lapse of a massive progenitor star involving a rebound shock 
in the stellar interior, expelling stellar materials and leav- 
ing behind a remnant compact object, such as white dwarfs, 
neutron stars, and black holes. 

This paper is arranged as follows. The general back- 
ground information in provided in Section 1. Section 2 gives 
an account for the hydrodynamic formulation of the model 
problem and previously known solutions. Section 3 treats 
the sonic critical curve. In Section 4, we present the novel 
'quasi-static' asymptotic solutions. We describe self-similar 
shocks and the relevant jump condition in Section 5. In Sec- 
tion 6, we show various semi-complete solutions either with 
or without a shock. Section 7 contains numerical examples 
for a rebound shock within a collapsing star. Finally in Sec- 
tion 8, we provide summary and conclusions. 



2 FORMULATION AND KNOWN SOLUTIONS 
2.1 Nonlinear Hydrodynamic Equations 



Our analysis and results are applicable to the vast region 
sufficiently far away from the central sphere surrounding 
various core activities. In the spherical polar coordinates 
(r,9,4>), the standard nonlinear polytropic hydrodynamic 
equations are in the forms of 



dp 1 d , 2 v 
+ — — [r pu) 



dt 
dM 



dr 







dM _ 

dt dr 

dM 2 
__ = 4vrr p , 
or 



( du 

p = up 1 



du 
dr 



dp _ GMp 
dr r 2 



(1) 



(2) 



(3) 



(4) 
(5) 



where u(r, t) is the radial bulk fluid velocity at radius r 
and time t, p(r,t) is the gas mass density, M(r,t) is the 
total mass enclosed within radius r at time t, G = 6.67 x 
10 -8 g _1 cm -3 s -2 is the gravitational constant, and p is 
the thermal gas pressure. The Poisson equation relating the 
gravitational potential $ (such that d^/dr — GM/r 2 ) and 
p is consistently satisfied. In equation of state @, the pro- 
portional coefficient /t is independent of t and is assumed 
to be a global constant. This n would be the same as the 
notation K(t) in e quation (11) of | Suto fc Si lk (1988) when 
the parameter n of ISuto fc Silkl Jl988) (and also in this pa- 
per to be defined presently) is equal to 2 — 7. The range 
1 < 7 < 4/3 for the polytropic index 7 is adopted in our 
current model analysis. 1 



2.2 A Self-Similar Transformation 

We introduce the self-similar transformation dSuto fc Silkl 
1988) with respect to a new dimensionless independent vari- 
able x as follows. 

r = ax , u = bv , p = ca , p = d(3 , M = em , (6) 
where the scaling factors a(t) through e(t) are defined by 



k 1/2 t" , 
kt 2 "' 4 

AtyG 



b = k 1 ' 2 ^- 1 



£.3/2^371-2 



AltGt 2 ' 



(3n-2)G ' (7) 

with k and n being two constants. Here v(x), a(x), @(x) and 
m(x) are dimensionless functions of x only. Substituting ex- 
pressions (|SJ and into the original polytropic hydrody- 
namic equations Q through JSJ and assuming n = 2 — 7 for 
a global constant k = fc(47rG) 7_1 , we obtain 



m — ax (nx — v) 

' 2 / 
a — a (n 



, , nx — v 
~ 1)U+ 3^2 Q 

2{x — v)(nx — v) 



\a(nx — v) 2 — 7a 7 ] 



(8) 



(9) 



and 



The present theoretical model problem treats the self- 
gravitational fluid dynamics under the spherical symmetry. 



1 See Appendix [X] for a brief discussion on cases of 7 ^ 4/3. 
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(n — l)av(nx — v) + 
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«) 2 2 

-a 



(3n - 2) 



27a 



\a(nx - 



■ya 



(10) 



So far , the development is the same as that of ISuto fc Silkl 
( 1988) with the requirement of n = 2 — 7. From equation 
0, it is clear that na; — v > is necessary for m > as a 
physical requirement for solutions. 

2.3 Known Similarity Solutions 

Various solutions of the two coupled nonlinear ordinary dif- 
ferential equations (ODEs) © and HUH have been known 
previously. These solutions with n = 2 — 7, which were first 
established in the isothermal case and then generalized to 
the polytropic case, are summarized below. 

The static solution of a singular polytropic sphere (SPS) 
is characterized by 



v = 0, 



2 7 (3n - 2) 



l/r, 



-2/r, 



2 7 (3n - 2) 



l/n 



3-2/n 
r 1 



(11) 



This is to be compared with the singular isothermal sphere 
(SIS; e.g., Shu 1977) and with the magnetized singular 
isothermal sphere (mSIS; Yu & Lou 2005). 

The Larson-Penston type of solutions is generalized to 

2x 2 
3 ' 



2(n-2/3) 3 
-x 



(12) 



3 3 
which was found by Larson (1969a, b) and Penston (1969a, 
b) in the isothermal case with n = 1; this is also known as 
the non-relativistic Einstein de-Sitter polytropic expansion 
solution (see Shu et al. 2002 for the isothermal counterpart) . 
The asymptotic solution finite at large x is given by 



nA 



Ax 



3n - 

-2/n 



+ 



27A 



7-1 



(n-2)/r. 



+ Bx 



(n-l)/r, 



(13) 



where A and B are two constants of integration. The 
isothermal count erpart solut ion with B = and n = 1 
was discussed bv IShul (jl9T7fl . and the more general solu- 
tion containing the free param eter B was first obtained by 
IWhitworth fc Summers! |l985) . In this paper, asymptotic so- 
lutions with B — and v > as the leading term at large x 
is referred to as breeze solutions, asymptotic solutions with 
B — and v < as the leading term at large x is referred to 
as the contraction solution, asymptotic solutions with B > 
are referred to as outflow or wind solutions and asymptotic 
solutions with B < are referrred to as inflow solutions. 

The leading diverging behaviour of central free-fall col- 
lapse solution in a polytropic gas at small x is characterized 
by 

1/2 

zm(U) 

v 



2m(0) 
(3n - 2)x 

(3n - 2)m(0) 
2x 3 



1/2 



(14) 



A Sample of Polytropic Solutions 
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Figure 1. Polytropic hydrodynamic similarity solutions with 
7 = 1.3 and n = 0.7 as examples of illustration for the counter- 
parts of the isothermal similarity solutions formerly constructed 
by various authors. The dash-dotted curve represents the sonic 
critical curve as defined in the text. The heavy dashed curve 
marked with 'Shul' is a solution with a contraction at large x and 
a free-fall collapse at small x without crossing the sonic critical 
curve. The heavy solid curve marked with 'Shu2' is the polytropic 
EWCS. The two heavy dotted curves labelled by 'WS1' and 'WS2' 
are inflow and outflow solutions at large x with free-fall collapse 
at small x, respectively. The light solid solution labelled 'LP' rep- 
resents a Larson-Penston type of polytropic solution. The heavy 
solid curve labelled by 'Hunter' is a Hunter type polytropic so- 
lution obtained by a phase matching procedure with the sonic 
critical point parameters (a, x, v) = (0.1641, 3.4865, 1.5711) 
and a central reduced density parameter a* = 22.038, and the 
heavy dashed curve labelled by 'LS' is a polytopic EECC so- 
lution first constructed by Lou & Shen (2004) for the isother- 
mal case with two sets of relevant sonic critical point parame- 
ters (a, x, v) = (164.69, 0.002978, -2.4497) and (a, x, v) = 
(0.1639, 3.4684, 1.5586). The last two curves happen to be fairly 
close to each other when x > 0.1. 



with m( x) approach ing a finite value m(0), which was first 
found bv lStnJ l|l977l ^ for the isothermal case with n = 1. 

The Hunter type of asymptotic solutions at small x are 
given by 



2 ai 1 ^ 
v = —x 

3 157 

„.( 2 -7) 



2 

Q* 

3 



67 



2 , 
x + ■ 



(15) 



where a* is an constant of integration. The is otherma l coun - 
terpart of this solution was first obtained bv lHunterl lll977f) 
with n = 1 and 7 = 1. 

These solutions are summarized here as the polytropic 
counterparts of the relevant results in the isothermal case. 
Figure Q is a collection of examples for previously known 
solutions with parameters 7 — 1.3 and n = 0.7. The EWCS 
constructed here is the polytropic c ounterpart of the isother- 
mal EWCS obtained bv lShul ll9_7J]) . Given the more general 
equation of state studied bv IChena il9TSf) . the isothermal 
EWCS represents only a subset of all possible EWCSs (see 
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Table 1. Relevant parameters for the known types of global sim- 
ilarity polytropic solutions constructed for n = 0.7 and 7 = 1.3 
(see Figure HI are summarized below. Parameters m(0) and a* 
are provided when the inner asymptotic behaviours are in the 
free-fall type or the Hunter type, respectively. 



type 



m(0) 



B 



Shul 


0.350 


0.5 





/ 


Shu2 


0.273 


0.5 


-1 


/ 


WS1 


0.374 


1 


1 


/ 


WS2 


0.700 


0.404413 





/ 


LS 


3.58xl0~ 3 


2.34 


3.78 


/ 


Hunter 


/ 


2.37 


3.82 


22.038 



also McLaughlin 
mation) . 



Pudritz 1997 in the logotropic approxi- 



3 THE SONIC CRITICAL CURVE 

The points where the denominators on the right-hand sides 
(RHSs) of both ODEs © and itTUIl vanish constitute the 
singular surface. On this surface, the travel speed of distur- 
bances relative to the local flow speed is equal to the local 
sound speed. Global solutions of the two ODEs and 1101 
cannot cross this singular surface smoothly, unless the in- 
tersection happens to (i) lie on the so-called sonic critical 
curve, on which both the numerators and denominators of 
the RHSs in the two ODEs vanish, and (ii) have the first 
derivatives of v and a with respect to x satisfying the crit- 
ical conditions at this intersection point (see Whitworth & 
Summers 1985 for the isothermal case). There are possibil- 
ities to go across the sonic critical curve with weak discon- 
tinuities (e.g., Whitworth & Summers 1985) or with shocks 
(Tsai & Hsu 1985; Shu et al. 2002; Bian & Lou 2005; Yu, 
Lou, Bian & Wu 2006; Lou & Gao 2006). The sonic critical 
curve and the critical conditions are derived below (see also 
McLaughlin & Pudritz 1997). 



3.1 Determination of the Sonic Critical Curve 

For an isothermal gas, the sonic critical curve can be ex- 
pressed in the form of v and a as functions of x JShdll977t 
lLou fc Shenll2004l) . For a polytropic gas, the sonic critical 
curve can be determined numerically for x and v values from 
specified a values. As only two of the three vanishing equa- 
tions of numerators and denominators are independent, we 
just need to consider the following two equations. 



(16) 



(17) 



/ *i2 7-1 

(nx — v) — 7a 
and 

, , - (nx — v) 2(x — v)(nx — v) 
( n -^ + a W^)- A -x L =°- 

For nx — v > and thus m > 0, we then obtain 

V = nx - y7a (7 ~ 1)/2 (18) 

and 



n - 1 + 



3n - 2 



v^ 1 



„(7-l)/2 _ 2 7 7-1 



n(n- l)x (19) 



to determine the sonic critical curve. Given a specific a value 
(note that a > by definition), the corresponding x values 
can be solved from equation 1)19^ as 



n - 1 + 



3n - 2 



n- 1 



3n 



8n(l - n) 



[2n(l 



(7-l)/2 



(20) 



We then obtain the corresponding v from equation l)18|l once 
q and x are known. 



3.2 Completeness and Asymptotic Behaviours 
of the Sonic Critical Curve 

According to equation 1191 . in order to achieve a positive x 
value, a must be larger than a critical value a c defined by 

1/2 



(3n-2){2[2n(l-n)] + 1 - n) 



(21) 

When q > a c , equation 1191 gives two different positive x 
roots. It is also routine to verify that when taking the '+' 
and ' — ' signs in equation 1201 . x monotonically increases 
and decreases with increasing a, respectively. 

The asymptotic behaviours of the sonic critical curve 
can be inferred from equation 1201 . When a approaches +00 
and taking the '+' sign in equation 1201 . we have 



(7+l)/2 



/[n(l-n)(3n-2)] 



(22) 



while taking the ' — ' sign in equation 1201 . we have 
x ~ 2^{3n - 2)q (7 " 3)/2 , 



(7-l)/2 



-V7 



-, -(7-1)7(3-7) 



2^7(3n - 2) 



.(23) 



Equations I22H and 1231 1 show the asymptotic trends of v and 
x versus a, and is valuable in determining the entire sonic 
critical curve. Our analysis here is c omplementary to th e 
analysis on the sonic critical points bv lSuto fc Silkl il988). 

Fig. |2]shows the sonic critical curves for different values 
of 7 and n with n = 2 — 7. This figure shows clearly the 
qualitative difference in the asymptotic behaviours between 
the isothermal and polytropic cases. 



3.3 Eigensolutions across the Sonic Critical Curve 

As the sonic critical curve is a singular curve for the two 
nonlinear ODEs @ and 1101 , it is essential to d etermine the 
solutio n behaviours in the vicinity of this curve. ISuto fc Silkl 
( 1988) analyzed the solution behaviours in the vicinity of the 
sonic critical curve and obtain 



(1 + 7 )tj' 2 + 



+2(27- 



n+ 1 -47 + 4(7 - 1)- 

x 



+ 2 



(3n - 2) 



2n - 47 + 4 



+ (fe§r +2(n+2 ^- 2) =° 



(24) 



and 
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Sonic Critical Curves with Different 7 Values 



si 
\ 

, % 7 = 1.2, n = O.S 
1 v, V>. ' 



-1.5 



-X 

7 = 1.3, re = 0.7 



7 = 1.2,re = 0.9 s v 



1.5 



2.5 



3.5 



Figure 2. The sonic critical curves for different pairs of 7 and 
n values. Here, n = 2 — 7 is adopted for a conventional poly- 
tropic gas. The 7 = n = 1 line (the dashed line) is a straight line 
v = x — 1 for the isothermal case (Shu 1977; Lou & Shen 2004). 
The major differences between isothermal and polytropic asymp- 
totic behaviours are the divergent behaviour at small x and the 
different slope of — v versus x at large x of the polytropic cases. 



- 2a(x — v) /x 
(nx — v) 



(25) 



where the data set of (a;, v, a) is along the sonic critical 
curve and v' and a 1 axe the first-order derivatives of v and 
a with respect to x across the sonic critical curve. In gen- 
eral, there exist two eigensolutions, each corresponding to 
a pair of first-order de rivatives v 1 and a', that can cross 
the sonic critical curve. IWhitworth fc Summersl (I1985T) con- 
structed solutions crossing the critical curve with first-order 
derivatives satisfying equations 1241 and 1251 for the isother- 
mal cas e , but with discontinuities in higher order derivatives. 
iHunterl il986l) promptly pointed out that such solutions in- 
volve weak discontinuities and may not be physically valid. 
In terms of global solutions, we shall mainly consider ana- 
lytically smooth solutions that cross the sonic critical curve. 



4 NEW 'QUASI-STATIC ASYMPTOTIC 
SOLUTIONS 

We now describe the new 'quasi-static' asymptotic solutions 
in the present model framework. 

4.1 Singular Polytropic Sphere (SPS) and 
Quasi-Static Asymptotic Solutions 

The static solution of a singular polytropic sphere (SPS) is 
well known and is given by 

-l/n 

(26) 







2 7 (3n - 2) 



-2/r, 



as a thermal-gravitational e quilibrium (Suto & Silk 1988; 
note that n = 2 — 7 here). EE] l|l977ll used this solution 
to construct the so-called expansion wave collapse solution 
(EWCS) for the isothermal case with n — 1. We realize 



that this solution can in fact serve as an asymptotic 'quasi- 
static' solution when x becomes sufficiently small. In other 
words, the leading terms of v(x) and a(x) are described by 
equation 1261 . yet there may exist higher orders terms to 
form an asymptotic series solution. 

To be consistent, we assume up to the second orders 



Lx* 



o(x) + ■ 



(27) 



-l/r, 



-2/n 



+ Aa + - 



2 7 (3n - 2) 

where L and K above, and N below are three constants, 



Act = Nx 



K-l-2/r, 



= o{a) 



(28) 

with the notation o(a) indicating Aa/a — > as x — > + . 
Here, K may be complex in general (see the analysis below) 
and Re(K) > 1 is required such that the additional terms 
in equations 127^ and I28jl are indeed higher order terms 2 as 
compared to solution 126H . 

Substituting these expressions into equations @ and 
HUH for small x, we obtain respectively 

n 2 (^ + l)iV r 2 



n{K-\)N = 

leading to 
K 2 - (4/r 
and 

TV 



2 7 (3n - 2) 



K + 2- - 

n 



1 -V r ' 



2 7 (3n - 2) 



3)K + 2 



n 2 (l + K) 



2 7 (3n - 2) 



-l/n 



(29) 



(30) 



(31) 



Equation I31H shows that once a K root is determined, the 
ratio N/L is then obtained. Of course, L and TV can also be 
both zero, corresponding to the static SPS solution 1261 . 
The two K roots of quadratic equation 13UH are 



K 



n 2 2 



24 
n 



16 

i) 2 



1/2 



(32) 



For Re(K) > 1, n < 0.8 is required. When n s£ 12- 8^2, K 
is a real number, while for 12 — 8\/2 < n < 0.8, K becomes 
a complex number. 



4.2 Quasi-Static Solutions of the 
First Kind with a Real K 

The 'quasi-static' solution with a real K is apparent and is 
confined to a narrow range of n values, i.e., n ^ 12 — 8y/2. 
There exist two sets of such solutions corresponding to the 
two K roots of equation I32H being both larger than 1. In 
such cases, the two asymptotic solution are of the form 



Lx 



2 7 (3n - 2) 



-l/n. 



-2/n 



2 Note that in solution 1261 . the first-order term of v is viewed 
to be Ox as a linear function of x. 
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+ 



2L 
n 2 (l + K) 



2 7 (3n - 2) 



-l/n 



K-l-2/r, 



(33) 



where the two K roots are given by expression 1321 and L 
is an arbitrary parameter. 

4.3 Quasi-Static Solutions of the 

Second Kind with a Complex K 

The 'quasi-static' solution with a complex K (i.e., when 12 — 
8^2 < n < 0.8) appears special. For a complex K = K\ + 
iKi with K\ and K 2 both being real, we should take the 
real parts of both sides of equations 12711 and 12811 for real 
v(x) and a(x). Then the v(x) solution appears as 

v = Re(Lx K ) = Re\ [ Lx I<1 exp(iK2 In a;)] 

= x Kl \Li cos(-fsT2 In x) — L2 sin(.K2 In s)] , (34) 

where L\ and L2 are real and imaginary parts of parameter 
L, respectively. By equation 1321 . we have K\ — 2/n — 3/2 
and K2 = (—1/4 + 6/n — 4/n 2 ) 1 / 2 respectively. The corre- 
sponding second-order term Aa is 



Aq = — 

77/ 



2 7 (3n - 2) 



-l/n 



l [Li(l + Kx) + L 2 K 2 ] cos{K 2 lnx) 

- [L 2 (l + Ki) - UK2] sm(K 2 lnz) j 



(35) 



Since the relative phase factor in both expressions 1341 and 
1351 can be adjusted by choosing the two free parameters 
L\ and L2, there is no loss of generality to adopt a positive 
value of K 2 . 



5 SELF-SIMILAR SHOCKS AND 

POLYTROPIC JUMP CONDITIONS 

Shocks form when faster flows catch up slower flows. In our 
model, it is of considerable interest to discuss self-similar 
shocks with the shock front 'fixed' in the self-similar profile. 
Across such a shock front, the upstream and downstream 
regions experience a change in the specific entropy leading 
to a change in k parameter in similarity transformation (JSJ 
and Q. 



5.1 Jump Conditions for Polytropic Shocks 

We denote the upstream (i.e., a fluid flows from this side into 
a shock front) physical variables with subscript 1, and the 
downstream (i.e., a fluid flows on this side away from a shock 
front) physical variables with subscript 2. The location of a 
shock front is denoted with a subscript s 



t x s 



(36) 



where k is for the upstream region, i.e. fci = k or K\ = k. For 
k 2 = K.X 2 or k2 — k\ 2 with A 2 being a scaling parameter, 
the similarity transformation in the downstream side is then 



r 2 = a 2 x 2 , 



U2 



&2«2 



pi 



C2012 



p 2 = (I2P2 



M 2 



e 2 m 2 , 



(37) 



where the five time-dependent scaling functions 0,2 to e 2 are 
defined explicitly by 



a 2 



d 2 = X 



• 1/2 t n 
kt 2n ~ 4 



xk 1 ' 2 ^- 1 



1 



C2 



e 2 = y 



^3/2^3n-2 



4ttG ' " (3n-2)G ' 

At the shock front, we should have n = T2 

X\ — X s = \X2 ■ 



(38) 

r s and thus 
(39) 



Given xi, Vi and ai, we need only to calculate the corre- 
sponding x 2 , V2 and 012 across a shock, and the parameter 
A is automatically determined by equation 1391 and hence 
the downstream transformation 1371 and 1381 is known. 

We now describe the jump conditions for shocks in a 
polytropic gas. Using the standard procedure, we choose the 
framework of reference in which the shock front is instantly 
at rest. The jump conditions in this shock framework of 
reference are then given below. The mass conservation is 



[p(u s - u)] * = 



(40) 



where u s = dr s /di is the radially outward travel speed of 
the shock front and, following the notational convention, we 
denote the difference of the corresponding physical variables 
between the upstream and downstream sides by enclosing 
them within a pair of square brackets with a superscript 2 
and a subscript 1. The momentum conservation is 

n2]2 



[rep 7 + p(u s - u) 2 ] = . 
The energy conservation becomes 



p{Us 



+ 



7k(w s 



(7- 



: 



(41) 



(42) 



Combining equations 1401 . 1411 and 1421 together with 
the self-similar transformation, we derive three shock condi- 
tions in terms of the self-similar variables 



ai(nx\ — vi) = Xa2(nx2 — V2) 



aj + ct\(nx\ — vi) 2 
2 7 



A 2 \a~l + a.2(nx2 — V2) 2 



(nil 



V!) 2 + 



(7- 



(nx 2 



V2Y + 



27 



(7-1) 



„7-l 



(43) 
(44) 



(45) 



5.2 Shock Solutions by the Jump Conditions 

By introducing new variables Yi = n — Xi/vi with i = 1, 2, 
we reduce equations 1431 through 1451 to 



(46) 
(47) 

(48) 



aiYi 


= 02T2 , 


^4 
x\ 


ail 1 = — 

x 2 


r 2 + 


2 7 aj' 1 


(7-1) A 



2 7 a]" 1 



J (7-1) x\ 
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New Type II Solutions Without Shocks 
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Figure 3. Semi-complete similarity solutions approaching type 
II 'quasi-static' asymptotic solutions at small x without shocks 
for 7 = 1.3 and n = 0.7. Here, — v versus x is shown. The dash- 
dotted line is the sonic critical curve. The solid solution is an 
inflow solution which is obtained by integrating from (x, v, a) = 
(0.236, —0.975, 1) on the critical curve inward; and the dashed 
curve is an outflow solution which is obtained by integrating from 
(x, v, a) = (1.81, 0.395, 0.17) on the critical curve inward. 
The two dotted straight lines are the abscissa and ordinate axes, 
respectively. For a sufficiently small x, both solutions approach 
the Type II 'quasi-static' asymptotic solution. 



New Type I Solutions Without Shocks 
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Figure 4. Semi-complete similarity solutions approaching type 
I 'quasi-static' asymptotic solutions at small x without shocks 
for 7 = 1.32 and n = 0.68. Here, — v versus x is shown. The 
dash-dotted line is the sonic critical curve. The solid curve is an 
inflow solution which is obtained by integrating from (x, v, a) = 
(0.0932, —1.086, 1) on the sonic critical curve inward; and the 
dashed curve is an outflow solution obtained by integrating from 
(x, v, a) = (1.52, 0.280, 0.07) on the sonic critical curve inward. 
The two dotted straight lines are the abscissa and ordinate axes, 
respectively. For a sufficiently small x, both solutions approach 
the Type I 'quasi-static' asymptotic solution. 



Figure 5. Semi-complete similarity shock solutions containing 
Type II 'quasi-static' asymptotic solutions for 7 = 1.3 and n = 0.7 
are shown here in terms of — v versus x. The heavy dotted curve is 
a shock solution with an inflow in the outer portion, constructed 
by integrating from (x, V, a) = (0.751, —0.426, 0.3) on the sonic 
critical curve with x S2 = 0.5; the heavy dash-dotted line is a 
solution with an outflow in the outer portion, constructed by in- 
tegrating from (x, v, a) = (2.47, 0.866, 0.16) on the sonic critical 
curve with x S2 = 2; the light dashed curve is a solution with a 
contracting outer portion, constructed with parameters A = 0.5, 
B = and a chosen x B1 = 1.5; the heavy solid curve is a solu- 
tion with a static outer portion (i.e., part of a SPS), constructed 



with parameters A = 0.404, B = and a chosen x B 



1.6; the 



heavy dashed curve is a shock solution with a breeze in the outer 
portion, constructed with parameters A = 0.402 and B = and 
choosing x ai = 1.5. The two light dotted straight lines are the 
abscissa and ordinate axes, respectively. For small x, all solutions 
approach the Type II 'quasi-static' asymptotic solution. 



We intend to solve for T2, Q2 and xi with known values 
of Ti, «i and x\. By substituting C12 = aiTi/T^ of equa- 
tion 1461 . solving for x\ from equations I47H and 1481 re- 
spectively, and eliminating x\ accordingly, we arrive at the 
following quadratic equation in terms of T2 



(7 + 1) ^2 

— o 1 2 

27 



K/zf+air?). 



+ 



air a 
(7-1) 

27 



r? + j7^7 



7-1 



(49) 



Excluding the trivial solution F2 
of equation 1491 in the form of 



Ti, we obtain one root 



r 2 = 



2 7 a7 



7-1 



+ 



(7- 



(7+1)^! ' (7 + I) 
In terms of the upstream Mach number M\ defined by 
(u s — Mi) 2 _ pi(u B — ui) 2 _ x\Y\ 



Ml 



expression 15U1 reads 

U2 I j 2 



7Pi 



7-1 ' 
7«i 



(50) 



(51) 



£2 



+ 



(7-1) 



Ui li (7+lVW 2 (7+1) ' 
this is precisely equation (89.6) of lLandau fc Lifshit j il95' 



(52) 
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Figure 6. Semi-complete similarity shock solutions containing 
Type I 'quasi-static' asymptotic solutions for 7 = 1.32 and n = 
0.68 are shown in terms of —v versus x. The heavy dotted line 
is a solution with an inflow in the outer portion, constructed by 
integrating from (a;, V, a) = (0.812, —0.243, 0.1) on the sonic 
critical curve with a chosen x S2 = 0.7; the heavy dash-dotted line 
is a solution with an outflow in the outer portion and constructed 
by integrating from (x, v, a) = (1.52, 0.280, 0.07) on the sonic 
critical curve with a chosen x S2 = 1-2; the light dashed curve is a 
solution with a contracting envelope, constructed with parameters 
A = 0.15 and B = and a chosen x B1 = 1.5; the heavy solid 
curve is a solution with a static outer portion of SPS, constructed 
with parameters A = 0.114 and B = and a chosen x s , = 1.3; 
the heavy dashed curve is a solution with a breeze in the outer 
portion, constructed with parameters A = 0.1135 and B = and 
a chosen x B1 = 1.3. The two dotted straight lines are the abscissa 
and ordinate axes, respectively. All solutions approach the Type 
I 'quasi-static' asymptotic solution at small x. 



It is straightforward to prove that the following set of 
inequalities are compatible with each other (M2 is the down- 
stream Mach number) 



r 2 



< 1 



U2 



A > 1 , 



U s — Ul 

k-2 > ki 



< 1 



Mi > 1 



K2 > Kl 



Ml < 1 



(53) 



We only give a proof for the compatibility between A > 1 
and M\ > 1 below as an example of illustration. First note 
that when .Mi > 1, we have z = < 1. From equation 

1471 . we then solve for A 2 = x\/x2 as 



A 2 = z 1 '[ 7 M?(1 - z) + l] = 



2 7 (1 



(7 + 1)* -(7-1) 



+ Z 1 



The first derivative of A 2 with respect to z gives 



dA^ 
dz 



7 (7 2 -!)(*-!) V- 
[( 7 + l)*_( 7 _l)] ! 



< 



(54) 



(55) 



for 7 > 1. Since A — ► 1 as z 1, it follows that when z < 1 
(i.e., M\ > 1), A > 1 holds and vice versa. Acco r ding t o the 
Zemplen theorem [see, e.g.. lLandau fc Lifshitzl £[.959)], for 
entropy to increase across a shock going from the upstream 
to downstream side, we have physical shock solutions with 
u s — «2 < u s — ui . While this general conclusion only holds 



for weak shocks, it is true in our case from equivalent in- 
equalities (15311 . That is, the upstream and downstream re- 
gions should be supersonic and subsonic, respectively. 



5.3 Solutions and Physical Interpretation 

For the convenience of analysis and comprehension, we 
adopt the solution procedure and presentation in the fol- 
lowing manner. We do not distinguish between xi and X2 in 
the unshocked regions, even though the self-similar transfor- 
mations are different for the upstream and downtream re- 
gions. We use x si and x S2 for the shock location under the 
similarity transformations of the upstream and downstream 
sides, respectively. The ratio A = x si /x 32 is then fixed by 
shock jump conditions. In figure illustrations of shock solu- 
tions, we juxtapose curves for —vi versus x\ and for — V2 
versus X2 in the same figure, with the implied differences 
in downstream and upstream transformations, and link the 
two shock points (— v si , x 31 ) and (— v S2 , x 32 ) by a straight 
line, indicating that these two points represent the demar- 
cation of upstream and downstream quantities on the two 
sides of the shock front, respectively. A reader of these fig- 
ures may interpret the curves to the left of this straight line 
(a; < x S2 ) as inner downstream solutions, and the curves to 
the right of this straight line (x > x si ) as the correspond- 
ing outer upstream solutions. As several shock solutions are 
pact together within one figure, one should also note that 
for different shock solutions, this demarcation is also differ- 
ent because x s values are different. The advantage of such 
presentations is that the critical curve of both regions co- 
incide in the figure, i.e., the curve provided in these figures 
serves as the critical curve in both regions. Furthermore, one 
can obtain the x si and x S2 values directly from the figure 
to compute A values accordingly. In our figures, the sear- 
ings of physical quantities are different in the presentation 
for upstream and downstream regions, and the scalings are 
only for the reduced quantities. For example, one unit in the 
two regions of the figure represents the same difference of x 
values, but not the same difference of r values; the situation 
is similar in terms of the correspondence between v and u 
values. 

Once we obtain a shock solution as described above, the 
procedure of obtaining a physical solution of our hydrody- 
namic model is to assign a k value of either region, i.e., ki or 
k2, and obtain the k value for the other region by the simple 
relation = feiA 2 . With the separate self-similar transfor- 
mations in the upstream and downstream regions, we readily 
compute the corresponding physical quantities of the solu- 
tion. In the dimensional form, the difference in the scaling 
of physical quantities due to the similarity transformations 
disappears. For example, r Sl = r S2 now, and a unit in the 
two regions in the figures represents the same difference of 
r values. 



6 SEMI-COMPLETE SOLUTIONS WITH 
THE 'QUASI-STATIC ASYMPTOTIC 
BEHAVIOURS 

Semi-complete similarity solutions include those that do not 
intersect with the singular surface, that go across the sonic 
point analytically, and that cross the sonic critical curve 
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Table 2. Parameters A, B and L of asymptotic behaviours for 
solutions without shocks presented in Figures [3l and l4l 



Table 4. Parameters A, B and L for asymptotic behaviours and 



shock parameters x s 



and A for similarity shock solutions 



n 


type 


A 


B 


L 


0.7 


outflow 


0.673 


0.852 


l.ll + 0.735i 


0.7 


inflow 


0.0840 


-2.27 


-11.4 + 19.7? 


0.68 


outflow 


0.176 


0.573 


0.452 


0.68 


inflow 


0.0142 


-2.55 


-301 



with type I 'quasi-static' asymptotic behaviours shown in Fig. IB1 



Table 3. Parameters A, B and L of asymptotic behaviours and 
shock parameters x si , x S2 and A for global similarity shock solu- 
tions approaching type II 'quasi-static' asymptotic behaviours at 
small x as shown in Figure IB1 



type 


A 


B 


L 


x S2 




A 


inflow 


0.0765 


-0.498 


-1.78 


0.7 


0.7009 


1.0013 


outflow 


0.184 


0.412 


0.452 


1.2 


1.225 


1.021 


contract 


0.15 





-0.203 


1.43 


1.5 


1.046 


static 


0.114 





-0.0196 


1.298 


1.3 


1.0018 


breeze 


0.1135 





-0.0248 


1.298 


1.3 


1.0015 



•Quasi-Static' Solutions at Small x 



type 


.4 


B 


L 


X S2 




A 


inflow 


0.227 


-1.15 


2.09 - 0.841i 


0.5 


0.511 


1.022 


outflow 


1.19 


1.56 


3.57+1.87i 


2 


2.07 


1.034 


contraction 


0.5 





0.289+0.188i 


1.46 


1.5 


1.024 


SPS 


0.404 





0.00233 - 0.0973i 


1.59 


1.6 


1.0058 v j 


breeze 


0.402 





-0.00207 - 0.0562i 


1.497 


1.5 


1.0018 



with shocks in the region < x < +oo (see Lou & Shen 
2004; Bian & Lou 2005; Yu et al. 2006). Either with or with- 
out shocks, semi-complete solutions approaching the 'quasi- 
static' asymptotic behaviours at small x are constructed in 
the present model. Specifically, we are able to construct self- 
similar shock solutions approaching the 'quasi-static' asymp- 
totic behaviours at small x yet with inflow, outflow, SPS, 
breeze or contraction in the outer envelope. In this section, 
we choose n = 0.7 and n = 0.68 for type II and type I 'quasi- 
static' asymptotic behaviours at small x, respectively. 



"2 5_ 



- ^ type II ( 


n =0.7} 








/ \ 1 
1 \ / ' 
/ \ / ' 
/ \ / ' 
J \l ' 








type I ( 


n =0.68) 





vj X 1 



6.1 Smooth Similarity Solutions without Shocks 

We are able to construct semi-complete similarity solutions 
without shocks simply by integrating inward from a certain 
point on the sonic critical curve, using the analytical eigen- 
solution with a smaller v' in the close vicinity of this cho- 
sen point [equations 12 H and 125H determine the eigensolu- 
tion]. Figures |3] and 2] are examples of illustration for semi- 
complete similarity solutions without shocks, approaching 
type II and type I 'quasi-static' solutions at small x, respec- 
tively. In Figure 01 only similarity solutions with smaller K 
values in equation 1321 are constructed by integrating from 
the sonic critical curve using one of the eigensolutions. 



6.2 Similarity Solutions with Shocks 

Based on the jump conditions, we are able to construct sim- 
ilarity shock solutions with inflow, outflow, breeze, SPS or 
contraction in the outer portion, connected to either type II 
or type I 'quasi-static' asymptotic solutions at small x. In 
constructing semi-complete shock solutions with inflow or 
outflow for the outer asymptotic behaviours, we use an in- 
ner solution that can cross the sonic critical curve smoothly 
by integrating inward from the sonic critical curve using an 
eigensolution, choose an x a2 value in the downstream region 
along this inner solution, determine the upstream quanti- 
ties by the shock jump conditions, and integrate outward 



Figure 7. 'Quasi-static' asymptotic solutions for small x val- 
ues. The presentation is in the v/x K versus In a; form for type I 
'quasi-static' asymptotic solution, and in the v/x Kl versus lnx 
form for type II 'quasi-static' asymptotic solution. The type I 
solution is obtained by integrating inward from (x, v, a) = 
(0.812, —0.243, 0.1) on the sonic critical curve, and the type 
II solution is obtained by integrating inward from (x, v, a) = 
(1.81, 0.395, 0.17) on the sonic critical curve. 



hand, in order to construct similarity solutions with breeze, 
SPS or contraction in the outer portions, we integrate in- 
ward from a large x value (say 10 5 ) to reach a certain value 
x 31 , determine the downstream quantities and integrate fur- 
ther inward to small x values. Only when the inner solutions 
approach the 'quasi-static' asymptotic solution do we suc- 
ceed in constructing such global shock solutions. Figures [K] 
and Hj] are examples of illustration for self-similar shock so- 
lutions approaching type II and type I 'quasi-static' asymp- 
totic behaviours at small x, respectively. In both Figs.|H|and 
Hj] examples of inflow, outflow, contraction, breeze and SPS 
for outer portions are illustrated. For type I 'quasi-static' 
asymptotic solutions, we again obtain solutions only with 
smaller K values in equation (13211 . 

6.3 Asymptotic Behaviours of Sample Solutions 

All solutions presented in this section approach the 'quasi- 
static' (either type I or type II) asymptotic behaviours at 
small x values. Figure [7| displays two examples illustrating 
such asymptotic behaviours in terms of v/x K versus In a; for 
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Figure 8. Sample shock solutions for polytropic gas dynam- 
ics. Solutions 1 and 2 are obtained by integrating inward from 
(x, v, a) = (2.533, 0.9074, 0.1598) on the sonic critical curve 
and choosing x B2 = 1.6198 and i, 2 = 1.64, respectively. Solu- 
tions 3 and 4 are obtained by integrating inward from (x, v, a) = 
(1.793, 0.4374, 0.067) on the sonic critical curve and choosing 
x S2 = 1.1 and x S2 = 1.2, respectively. The relevant parameters 
for asymptotic behaviours and for shock conditions are summa- 
rized in Table The two dotted straight lines are abscissa and 
ordinate axes, respectively. The straight lines in the middle rep- 
resent shocks, and the two x B values correspond to the same r s . 



Table 5. Parameters for asymptotic solutions and shocks in the 
polytropic shock solutions displayed in Figure 151 



No. 


1 


2 


3 


4 


7 


1.3 


1.3 


1.32 


1.32 


A 


3920 


23.5 


5.40 


0.493 


B 


-0.00442 


0.562 


-0.792 


0.371 


L 


3.83 + 2.27i 


3.84 + 2. 28i 


0.441 


0.441 


x si 


29.79 


5.013 


3.288 


1.578 


X S2 


1.6198 


1.64 


1.1 


1.2 


A 


18.39 


3.056 


2.989 


1.315 



type I solution and v/x Kl versus In a; for type II solution. For 
the type II 'quasi-static' solution, the sine (cosine) pattern 
is apparent, while for the type I 'quasi- static' solution, the 
tendency to a constant value is also fairly clear. For the 
convenience of comparison and reference, we also list all the 
relevant parameters x 31 , x S2 and A for shocks and A, B and 
L for asymptotic behaviours in Tables [5] [3] and [I] 



core collapse and a rebound shock scenario for a supernova 
explosion. 

At the onset of gravitational core collapse inside a mas- 
sive star, the central region suddenly loses pressure and 
a nearly free infall towards the center ensues; meanwhile, 
the information of core collapse travels outward through 
the stellar interior to reach the envelope and a self-similar 
collapse phase grossly characterized by asymptotic free-fall 
solution 114H at small x oc r/t n may gradually emerge. 
By physical estimates, this happens on a timescale of a 
few seconds. For small x, we obtain 3 the radial flow speed 

U OC _ t (3«/2)-l r -l/2 i the mass density p x t (3n/2)-2 r -3/2 j 

the enclosed mass M fc 3/2 m(0)t 3n_2 /[(3n - 2)G], and 
the core mass accretion rate M — > fc 3/2 m(0)t 3(n_1) /G with 
n > 2/3. It should be noted that before the onset of core 
collapse, the progenitor star is expected to have a wind in 
general and the star itself can engage in global oscillations, 
be they acoustic p— modes or internal gravity g— modes. For 
example, for purely radial stellar acoustic oscillations, the 
stellar envelope and interior can move radially either out- 
ward or inward depending on the phase of oscillation during 
the rapid core collapse. These plausible physical conditions 
may lead to different self-similar evolution behaviours during 
a free-fall core collapse. Also, the material degeneracy can 
already set in the core well before the onset of gravitational 
core collapse. Moreover, the mass density divergence asso- 
ciated with the free-fall collapse solution 11411 can rapidly 
trigger core degeneracy. In reality, no inifinity should arise. 

As mass accumulates towards the center, singularity 
arises mathematically and similarity disappears. As the de- 
generate core pressure rapidly builds up to resist the ma- 
terial infall, a powerful rebound shock emerges around the 
core and drives out most gas materials as it ploughs through 
the infalling stellar envelope, leaving behind a remnant com- 
pact object (such as a neutron star or a black hole 4 ) with 
high mass density in the degenerate core. For instance, for 
the formation of a neutron star, the breakout of such a re- 
bound shock through the stellar photosphere is crucial to 
initiate a supernova explosion observed at terrestrial ob- 
servatories (e.g., Bethe 1995; Burrows 2000 and Lattimer 
& Prakash 2004 for reviews of supernova and neutron star 
physics and extensive references therein). Again, we pro- 
pose that after the emergence of a rebound shock around 
the collapsing core, a self-similar evolution grossly charac- 
terized by either 'quasi-static' asymptotic solution 13311 or 
quasi-static' asymptotic solution 1341 and 1351 may gradu- 
ally appear and persist as the rebound shock travel outward 
against infalling materials. The self-similar evolution phase 
of a rebound shock may last on timescales in the range of 
~ 10 4 - 10 5 s or longer. 



7 SELF-SIMILAR REBOUND SHOCKS 
IN CORE COLLAPSING STARS 

7.1 Physical Scenario 

The core collapse of a massive star involves a multitude of 
physical processes of nuclear burning, radiation pressure, de- 
generate materials, electrons, and neutrinos (Hirata et al. 
1987; Bionta et al. 1987; Bethe 1993). Here, we only con- 
sider several gross aspects from the point of view of fluid 
dynamics involving self-gravity and try to outline a free-fall 



3 Bethe (1993) indicated u oc r -1 / 2 and p oc t~ 1 r~ 3 / 2 outside 
the (rebound) shock and after the start of (core) collapse. This 
would correspond to an index n = 2/3; for a conventional poly- 
tropic gas with n = 2 — 7, this would in turn imply 7 = 4/3 for an 
extremely relativistically hot gas (e.g., Goldreich & Weber 1980). 
We take n > 2/3 in our model framework to make the enclosed 
mass M physically meaningful in self-similar transformation 171 . 

4 We will also discuss the possibility of producing or exposing 
a white dwarf in a similar process separately. For a low-mass 
progenitor star, the formation of a proto-white dwarf in the core 
may or may not involve a violent rebound shock explosion. 
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For 'quasi-static' asymptotic solution 133H of the first 
type at small x, we have u = Lk^~ K)/ H~ nl ^ K - 1) ~ 1 r K , 

— 2/n | ,—n(K-l) K—l — 2/n 11 j 3—2/n j ■> 

pocr'+rf 'r , M oc r , and M oc 

Lt _ n(if _ 1) _ l7 , K+2 _ 2/n ^ where L < Q for an infloW) K > Xj 

and e is a sufficiently small parameter. In other words, the 
radial flow speed u decreases with time and approaches zero 
as r — * 0; the mass density p diverges as r — > to induce 
material degeneracy; being independent of time t, the en- 
closed mass M vanishes as r — > 0; and the core mass ac- 
cretion rate M decreases with time and approaches zero as 
?- — > by quadratic equation (13211 . A very important and 
testable prediction of this first type of self-similar rebound 
shock explosions is that the mass density p for stellar ma- 
terials scales as r~ 2 ^ n with 2/n < 3 behind the rebound 
shock after a sufficiently long lapse of time; for a conven- 
tional polytropic gas with n — 2 — 7 and 7 ^ I, we would 
also have 2 < 2/n. Ideally, this mass density profile can be 
better preserved with a more or less quiet central nonro- 
tating neutron star or Schwarzschild black hole, otherwise 
central activities or winds inevitably destroy central part of 
the density profile. An immediate example comes to mind is 
the Cassiopeia A supernova remnant which appears more or 
less spherical with a nonrotating X-ray bright neutron star 
left behind (i.e., no X-ray pulsations are detectable so far by 
Chandra satellite observations). 

For 'quasi-static' asymptotic solution I34H and 1351 of 
the second type (with logarithmic oscillations) at small 
x = r/(k 1 ^ 2 t n ), we then have the radial flow speed u = 
jfe( i_jf 1 )/2 t -„(if 1 -i)-i r .Jf 1 [ ii cos^lnx) -L 2 sm{K 2 \nx)], 
the mass density p oc r~ 2/n + ni (5n - 4)/2 r-- 5/2 {[Li (1 + 
A'i) + L 2 A" 2 ] cos(A 2 lnx)-[L 2 (l + K 1 )-L 1 K 2 ] sin(A 2 In a;)}, 
the enclosed mass M oc r 3 - 2/n , the core mass accre- 
tion rate M oc fc( 1 -- K i)/ 2 t-"( Jc i- 1 )- 1 r 1/2 [L 1 cos(K 2 In a?) - 
L 2 sin(A 2 Ins)], where n is a readily identifiable coefficient, 
Li and L 2 are two small parameters, K\ — 2/n — 3/2, 
K 2 = (-1/4 + 6/n - 4/n 2 ) 1/2 and n > 2/3. In reference 
to the first type of self-similar 'quasi-static' asymptotic so- 
lutions 1331 , the qualitative distinction of this second type of 
self-similar 'quasi-static' asymptotic solutions 1341 and 1351 
is of course the oscillatory feature. The overall strength of 
this oscillatory feature tends to die out with increasing time 
t. Nevertheless, at a given time t, the mass density profile of 
p may still retain a considerable oscillatory feature especially 
for small r or small x. Again, we propose to use the data of 
the Cassiopeia A supernova remnant or other suitable su- 
pernova remnants to identify or search for such oscillatory 
features. A series of our self-similar shock solutions with 
'quasi-static' asymptotic solutions of first and second types 
are consistent, at least qualitatively, to the major features 
of such a plausible rebound shock scenario for supernova 
explosions. 

To be more specific and tangible, we present a series 
of self-similar rebound shock solutions approaching central 
'quasi-static' asymptotic solutions in Figure |H] They serve 
as examples of illustration for a possible hydrodynamic sim- 
ilarity model of a rebound shock in a core collapsing pro- 
genitor star. Among various physical possibilities, these re- 
bound shocks are constructed as such, because we would 
like to demonstrate similarity shock solutions with (i) an 
inflow right on the upstream side of the shock, which ac- 
counts for the creation of a rebound shock propagating out- 
ward through a collapsing envelope; (ii) a self-similar shock, 



modelling a rebound shock after encountering a highly com- 
pressed over-dense core in a gravitational collapse; (iii) a 
massive outward expansion behind (downstream of) the re- 
bound shock, such that gas materials of the massive pro- 
genitor star are driven out by such a powerful bouncing; 
and (iv) a 'quasi-static' asymptotic behaviour at small x, so 
that eventually the dense inner core comes to a near equilib- 
rium as time goes on. These four features are qualitatively 
consistent with the rebound shock scenario of a perceived 
supernova explosion, leaving a compact object behind. We 
show shock solutions with both type I (solutions 3 and 4) 
and type II (solutions 1 and 2) 'quasi-static' asymptotic so- 
lutions at small x, because both satisfy the above intuitive 
physical requirements qualitatively. Solutions 1 and 3 have 
inflowing outer portions, while solutions 2 and 4 have out- 
flowing outer portions to encompass various conceivable sit- 
uations. In principle, there may also exist contracting and 
breezing outer portions with other features much alike the 
core collapse situations just described above. 

7.2 Spatial and Temporal Cutoffs 

For a massive progenitor star, the self-similar evolution of 
a core collapsing process with a rebound shock cannot per- 
sist all the time and thus a temporal or spatial cutoff needs 
to be specified. However, the similarity model is character- 
ized by a continuous mass density profile with no natural 
cutoffs present in solutions themselves. We thus refer to 
astronomical data to estimate a sensible boundary for our 
model application. In this model, we set a radius 7^ (for in- 
stance, r, = 10 6 cm if the compact object is a neutron star 
or black hole, and rt = 10 8 cm if it is a white dwarf) for our 
sample solutions; and we choose a radius r (for example, 
r = 10 12 cm) as the outer boundary for application of our 
sample solutions. A typical massive star has a radius of the 
order of ~ 10 12 cm (see Herrero et al. 1992 and Schonberner 
& Harmanec 1995 for observational data on masses and radii 
of main-sequence O and B stars). In principle, we could in- 
voke even larger stellar radii for even more massive stars. 

For these spatial cutoffs, the cutoff time is introduced 
below accordingly. A rebound shock emerges roughly around 
the core of a massive star when the collapse produces an 
over-dense degenerate core. It travels outward and evolves 
into a self-similar phase. We take the time when the shock 
crosses the inner radius (say, ~ 10 6 cm for a neutron star 
radius) as the initial time to apply our shock solution. This 
time would then correspond to 

Here, t 1 is obtained assuming our model to be valid from the 
beginning at t — 0. This estimate also roughly represents the 
time needed for the shock to evolve from the center to the 
inner reference radius r^. Specific estimates of fi are summa- 
rized in Table |7| with adopted parameters indicated. Since 
our 'quasi- static' asymptotic solution approaches a static 
state as t — -> +00, we do not set an upper bound in time for 
our rebound shock model. 

5 The progenitor star may oscillate radially prior to and during 
the onset of a gravitational core collapse as already noted earlier. 
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Figure 9. The negative radial flow speed profile — u at different 
time t values of the numerical example, with r in the logarithmic 
scale, for the model 1 in Figure Here, ti = 6.15 X 10 _5 s is the 
time when the shock crosses the inner radius rj, t2 = 2.29 X 10 4 s 
is the time when the shock crosses the outer boundary r Q ; t m \ = 
0.1s and t m 2 = 10 5 s are intermediate times between t\ and f-2, 
and between ti and t = oo, respectively. Negative radial speed 
profiles have been multiplied by various numerical factors shown 
in the figure for the compactness and clarity of presentation. 



7.3 Properties of Rebound Shock Solutions 

We note a common property of self-similar rebound shock 
solutions shown in Figure [g] As x — > +oo and as x — * + , 
the profiles of radial flow speed u and mass density p both 
become independent of time t. For a specified radial range 
of r within a star, the stellar configurations are either static 
or stationary both at the beginning (t — + + ) and at the end 
(t — ► +oo). In particular, the enclosed mass is given by 



M = 



nk 1/n A 



3-2/r, 



(3n- 2)G 

for asymptotic solutions at large x, whereas 

l/n 



M 



(3n - 2)G 



fc27(3n - 2) 



3-2/„ 



(57) 



(58) 



for the 'quasi-static' asymptotic solutions of both types at 
small x. 

These analytical results may be utilized to grossly esti- 
mate mass variations of a star during the core collapse and 
supernova explosion processes. As t approaches + , if one 
takes our model to be valid from the beginning, then the 
enclosed stellar mass M is known once a progenitor stel- 
lar radius r is specified. Likewise, if one takes our model 
to be valid after a sufficiently long time (i.e., t — + +oo), a 
quasi-static core equilibrium is ultimately achieved and the 
enclosed mass M is known for a specified core radius r^. 

A summary of the estimated enclosed mass values for 
both the outer stellar radius r and the inner reference ra- 
dius Ti for solutions presented above are given in Table HJ 
for the four rebound shock solutions in Fig. |H| We refer to 
the enclosed mass at t — > + as the initial mass (hence the 
subscript i n i), and the enclosed mass at t — > +oo as the ulti- 
mate mass (hence the subscript u u), respectively. The initial 



Table 6. The adopted k values and the initial and ultimate en- 
closed mass values for the four models shown in Figure [5] The 
subscript o denotes the outer enclosed mass, and subscript i de- 
notes the inner enclosed mass; while subscript i„i represents the 
initial mass, and „; t represents the ultimate mass. 



No. 


1 


2 


3 


4 


fci(cgs units) 


8.87 x 10 14 


3.21 x 10 16 


3.36 x 10 16 


1.735 x 10 


Ai2(cgs units) 


3 x 10 17 


3 x 10 17 


3 x 10 17 


3 x 10 17 


Mi,in«(M©) 


3.35 


3.38 


3.11 


3.18 


M iMt (M Q ) 


1.42 


1.42 


1.65 


1.65 


M o , ini (M ) 


24.1 


24.3 


7.02 


7.17 


M 0iuH (M e ) 


10.2 


10.2 


3.71 


3.71 



mass may roughly represent the mass of a progenitor star 
just evolved into a self-similar phase after the initiation of 
a core collapse, and the ultimate mass may roughly corre- 
spond to the mass of the remnant compact object as a long 
time has elapsed. The determination of the k values in both 
the upstream (subscript 1) and downstream (subscript 2) 
regions are described in the next subsection; while the inner 
and outer radii are set to be ri = 10 6 cm and r = 10 12 cm 
here for considering remnant neutron stars or black holes. 



7.4 Estimates for values of k Parameter 

Clearly, once the numerical value of the k parameter in sim- 
ilarity transformation JSJ and J7J is determined for a star, 
the shock model is then specified. Since 

k - 



pT(4 7 rG)T- 1 
Nk B T 



k B T 



Pt(4ttG)t- 



mpT- 1 (47rG)^ 



(59) 



where iV is the particle number density, m is the mean 
molecular (atomic) mass of the gas particle, and the sec- 
ond equality holds only for an ideal gas, we may in principle 
determine k parameter using the physical quantities in equa- 
tion I59H . Nevertheless, even for long-studied main-sequence 
stars, the determination of physical parameters in the stel- 
lar interior has been challenging, not to mention the relevant 
parameters for massive collapsing stars, which involve var- 
ious rapid physical processes and are extremely rare to be 
caught in action. Currently, some theoretical estimates of 
such parameters are available, although most of these pa- 
rameters are estimated for static configurations, e.g., main- 
sequence stars or neutron stars. These two stellar configura- 
tions as part of our model considerations, are the progenitor 
star and the remnant compact object left behind a rebound 
shock, respectively. While the value of k parameter is by no 
means constant during complicated processes before and af- 
ter the similarity evolution involving a constant k, we would 
presume for simplicity that the k value does not change sig- 
nifica ntly during the dynamical evolution. 

lAppenzeller fc Tscharnuterl Jl974l) 

iBond. Arnett fc Carrl (|l984l) iNadvozhin fc Razinkoval 
( 2005) and Sch aller et alJ il992T) gave theoretical estimates 



for the hydrogen-burning phase of massive stars with 
p c ~ lg/cm 3 and T c ~ 10 7_8 K. Using m = 0.5m p with 
m p being the proto n mass , we obtain approximately 
k ~ 10 17 cgs units. lArnettJ jl977f) . Bond et al. (1984), 
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Figure 10. The total enclosed mass at different time t values of 
the numerical example in logarithmic scales, for the first model 
in Figure IS1 Here, t± = 6.15 X 10 — 5 s is the time when the shock 
crosses the inner radius r<, ta = 2.29 X 10 4 s is the time as the 
shock crosses the outer boundary r a ; t m \ = 0.1s and t m 2 = 10 5 s 
are intermediate time values between t\ and f,2, and between t-2 
and t = oo, respectively. 

and ISchaller"et al. ( 1992) provided theoretical predictions 
for the late evolution phase of massive stars after the 
hydrogen-burning phase, and the results are p c ~ 10 8 g/cm 3 
and T c ~ 10 9 K. Using these late-phase parameters and 
m = 4m p /3 for helium atoms, we derive k ~ 10 16 cgs units. 
Also from the equation of state for degenerate neutrons, 
the corresponding p c ~ 10 35 dynes/cm 2 and p c ~ 10 15 g/cm 3 
may be appropriate for neutron stars (e.g., Shapiro & 
Teukolsky 1983), and the estimated k value is accordingly 
~ 10 17 cgs units. 

In the rebound shock model construction as we intend 
to discuss applications to SNe and the formation of remnant 
neutron stars, we may first assign k<i values for a reason- 
able enclosed mass range around ~ 1M§. To be consistent 
with estimates for typical neutron stars, the &2 value is es- 
timated to be ~ 3 x 10 17 cgs units for both type I and type 
II 'quasi-static' asymptotic solutions at small x. The value 
of fei is derived accordingly for each shock model, using the 
A parameter from the shock jump conditions. The k val- 
ues thus determined in our shock models are summarized in 
Table |S| We see that all the k\ values are of the orders of 
~ 10 15 ~ 17 cgs units, fairly close to the estimates for massive 
progenitor stars. 

7.5 Shock Model Presentation 

Figs. 101 and llOl as well as Figs. lll| and ll2| present the profiles 
of the negative radial flow speed —u and the enclosed mass 
M at different temporal epochs, respectively. 

Besides the timescale ti as the initial time of our model, 
we introduce another typical timescale ti which is the time 
when the shock reaches the outer radius r , namely 

f \ 1/n 



Figure 11. The radial speed profile at different time t values of 
the numerical example in logarithmic scales for the fourth model 
in Figure 181 Here, t\ = 7.18 X 10 — 5 s is the time as the shock 
crosses the inner radius r,, t2 = 4.78 X 10 4 s is the time as the 
shock crosses the outer boundary r D ; t m \ = Is and f m 2 = 10 6 s 
are intermediate times between t\ and *2, and between t2 and 
t = oo, respectively. The speed profiles have been multiplied by 
various factors as shown for the compactness and clarity of the 
presentation, 



Table 7. Below all time t in seconds and all enclosed mass in unit 
of the solar mass Mq. Here, t\ for the time as the shock reaches 
the inner boundary and t2 for the time as the shock reaches the 
outer boundary and the corresponding Mi and M for the four 
models in Figure GO Enclosed masses for £ m i and t m 2 as indi- 
cated in relevant figures are also included. The k parameters are 
according to Table IB1 



No. 


1 


2 


3 


4 


h 


6.15xl0~ 5 


6.04xl0~ 5 


8.16xl0~ 5 


7.18xl0~ 5 


t% 


2.29X10 4 


2.25 xlO 4 


5.44 xlO 4 


4.78X10 4 


Ma 


3.53 


3.54 


3.22 


3.23 


M a 


1.43 


1.43 


1.69 


1.69 


M ol 


24.1 


24.3 


7.03 


7.17 


M o2 


25.4 


25.4 


7.24 


7.28 


M iml 


1.06 


1.06 


1.85 


1.85 


M im2 


1.43 


1.43 


1.68 


1.68 


M oml 


24.1 


24.3 


7.03 


7.17 


M om 2 


15.6 


15.6 


5.78 


5.78 



We presume our model to be valid since t—*0. The specific 
t2 values relevant to various shock models are also contained 
in Table with adopted parameters clearly indicated. As a 
rebound shock cannot be seen from the Earth until the shock 
reaches the stellar photosphere and our outer boundary is 
set for the radius of a massive star, this time value roughly 
corresponds to the time for a rebound shock to travel out 
and be seen. 

Quantities at five different epochs: ti, £2, t = 00 and 
two intermediate times t m i and t m i are shown in the radial 
range n = 10 6 cm < r < r a = 10 12 cm in Figs. 151 111711111 
and 1121 The corresponding radial flow speeds of model 2 
are like those of model 4, and the corresponding radial flow 
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speeds of model 3 are like those of model 1. Meanwhile, the 
corresponding enclosed mass distributions of model 2 are 
like those of model 1, and the corresponding enclosed mass 
distributions of model 3 are like those of model 4. 



7.6 Analysis of the Rebound Shock Model 

We now describe model 1 results as an example of illustra- 
tion. Relevant data of this model as well as other models are 
summarized in Tables|^|and[2| Initially, the stellar interior is 
in a core collapse (inflow) stage, which is triggered by the ex- 
haustion of the central nuclear fuel. The initial profile is also 
significantly over-dense as compared to the eventual quasi- 
static configuration after a long lapse in time. The outer 
envelope remains collapsing as the rebound shock emerges 
around the central degenerate core and travels outward. The 
shock travels faster than the sound speed and the pertur- 
bation information cannot reach the outer region until the 
shock reaches there. At time ti, the enclosed mass within r 
is 24. 1M© and within rt, it is 3.53Mq. 

The rebound shock then breaks out from the stellar in- 
terior, travelling at a fantastic speed (mostly in the range 
of 10 7 ~ 10 9 cm/s). During the first few seconds, the core 
mass decreases rapidly. At time t m \ = 0.1s, the enclosed 
mass within n has reduced to I.O6M0, while within r the 
enclosed mass is still 24. 1M© because the time difference 
is small and the mass accretion is not significant. When 
the shock reaches the outer boundary, a slight accretion has 
changed the total enclosed mass within r to 25. 4Mq, yet 
the inner core has already experienced a minor accretion 
due to radial speed oscillation of our type II 'quasi-static' 
asymptotic solution, and the enclosed mass within the inner 
radius n has now increased to 1.43Mq. 

After the rebound shock passes over the outer refer- 
ence radius r , it is now the time for the massive star to 
drive out a significant amount of mass. As the drive grad- 
ually weakens, the total enclosed mass within r D at time 
tm.2 = 10 5 s falls to 15.6M0, while within Ti the total en- 
closed mass remains 1.43M0 as before. Eventually the en- 
closed mass within r becomes 10. 2Mq, and within Vi, it 
becomes 1.42Mq. From the core mass variation and the re- 
maining core mass within the inner reference radius r, in the 
end, we reasonably conjecture that the remnant core may 
become a proto-neutron star. The rebound shock velocity 
also conforms with usual estimates. 



7.7 Mass Ratio Evolution 

One interesting aspect of our rebound shock model is the 
ratio of the initial outer mass M ,mii roughly corresponding 
to the initial mass of a prescribed progenitor star, and the 
ultimate inner mass Mi :U u, roughly corresponding to the 
mass of a remnant compact object in the core. By equations 
15711 and 1581 . this mass ratio is simply 



Enclosed Mass Profiles at Different Epochs (4) 



M 0; i n i A 

Mi u it ~ A 3 ^ 



2 7 (3n - 2) 



1/n / s (3-2/n) 
To 



(61) 




In expression 1611 . the last factor involving the ratio r /ri 
which is somewhat arbitrary, although we can estimate 
proper radii for different progenitor stars and ultimate com- 
pact core objects in a sensibel range. The ratio of the two 



Figure 12. The total enclosed mass at different time t values of 
the numerical example, in logarithmic scales for the first model 
in Figure 151 Here, t\ = 7.18 X 10 — 5 s is the time when the shock 
crosses the inner reference radius r*, t% = 4.78 X 10 4 s is the time 
as the shock crosses the outer boundary r D ; t m \ = lis and t m 2 = 
10 6 s are intermediate time values between ti and <2, and between 
t2 and t = 00, respectively. 



radii naturally affects the mass ratio, but as n approaches 
2/3, or equivalently 7 approaches 4/3, this last factor would 
approach unity. The two factors together in front of this 
somewhat arbitrary last factor in expression 16H is the ra- 
tio of the initial mass and the final mass within the same 
radius (i.e., the last factor is 1 in this case). This can vary 
when 7 or equivalently n varies, but it also depends on x S2 
or (xo ,vo ,oto), which is the point on the sonic critical curve 
from which we construct the 'quasi-static' portion of our re- 
bound shock solutions (see subsection 6.2). Since the relation 
of A and A is explored numerically, we here only provide the 
corresponding results for several cases with n = 0.68 (see 
Table [SJ, and the cases with other n or 7 values are fairly 
similar. 

According to numerical results of Table |H] mass ratio 
<61> increases for smaller cto or for larger x 32 . This is con- 
ceivable, in the sense that as cto becomes smaller, the 'quasi- 
static' portion of the rebound shock solution has a higher 
v for more mass being driven out; and when x„ 2 becomes 
smaller, the outer portion has a higher v for more mass 
being driven out. Referring to Table |H] and the relevant def- 
initions therein, we also infer empirically that ratio 1 may 
approach 1, but whether this ratio 1 can be arbitrarily large 
is unclear. It appears that the choice of self-similar rebound 
shock solutions affects the mass ratio in a significant man- 
ner. We infer from this feature that in reality the mass ratio 
may also vary, and thus a compact object with a definite 
mass may possibly come from progenitor stars with differ- 
ent initial masses. 



8 REMNANT PROTO- WHITE DWARFS 

Conventionally, one does not invoke a rebound shock explo- 
sion to expose the degenerate core, a proto-white dwarf, of a 
progenitor star. It is a general belief that various mass loss 
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Table 8. Evolution of mass ratio M ,i n i/^i,ult f° r various pa- 
rameters when 7 = 1.32 and n = 0.68 are specified. For each value 
of ao, the point from which we construct our rebound shocks, we 
explore two cases, namely, rebound shock solutions with contract- 
ing and expanding envelopes: in the neighbouring two sets with 
the same ao values, the one with smaller x S2 gives an expanding 
envelope, and the one with smaller x B2 gives a contracting en- 
velope. In this Table, ratio 1 is the first two factors together in 

expression fiflV namely, A{n 2 /[\ 2 2-y(3n - 2)] } , and ratio 2 

__ / s(3-2/n) 

is the last factor in expression 1611 . namely, (r /r; ) 



ao 


X S2 


A 


A 


ratio 1 


ratio 2 


0.067 


1.2 


1.315 


0.493 


1.933 


4.357 


0.067 


1.1 


2.989 


5.401 


1.893 


4.266 


0.070 


1.1 


1.064 


0.2048 


1.497 


3.373 


0.070 


0.9 


1.737 


0.8338 


1.442 


3.250 


0.072 


1 


1.0628 


0.1808 


1.326 


2.989 


0.072 


0.8 


1.843 


0.879 


1.277 


2.878 


0.074 


1 


1.0248 


0.149 


1.216 


2.742 


0.074 


0.8 


1.262 


0.2662 


1.178 


2.655 


0.076 


1 


1.0100 


0.1326 


1.130 


2.546 


0.076 


0.8 


1.1238 


0.1766 


1.099 


2.478 


0.078 


1 


1.00383 


0.1219 


1.0575 


2.384 


0.078 


0.8 


1.0663 


0.1422 


1.0329 


2.328 


0.079 


1 


1.0000235 


0.1175 


1.0308 


2.323 


0.079 


0.8 


1.0496 


0.1320 


1.00440 


2.264 



mechanisms through winds are responsible for eventually 
producing or exposing proto-white dwarfs. However, given 
our theoretical framework outline in this paper, it is not ob- 
vious why a rebound shock cannot occur to blow away stellar 
envelope and produce a proto-white dwarf. After all, a pro- 
genitor star of mass ~ 6 — 8Mq has to throw away a mass 
amount of ~ 4.6 — Q.QMq in order to expose a proto-white 
dwarf in the core. 

To consider the formation of a proto-white dwarf in the 
degenerate core of a progenitor star in the high-mass end 
(say, ~ 6 — SMq), we would set r% to be a typical radius of 
several thousand kilometers (say, ~ 3 x 10 8 cm) for a white 
dwarf and r ~ 10 13 cm to be the stellar radius to enclose 
a mass in the range of ~ 6 — 8Mq for the corresponding 
progenitor star. After the initiation of core collapse and as 
a powerful self-similar rebound shock travels outward, the 
masses enclosed within r% and r decrease, sometimes oscil- 
late, and eventually approach finite values with the increase 
of time t. If the final mass M i>u i t within r, is less than the 
Chandrasekhar mass limit of 1.4M©, we would say that a 
proto-white dwarf is produced in the remnant core. Other- 
wise if the final mass Mj iU ( t within n is greater than 1.4M©, 
say ~ 2 — 3Mq , we would then hypothesize that dynamical 
instabilities associated with a proto-white dwarf exceeding 
the Chandrasekhar mass should have already happened be- 
fore reaching the final model state of t —* +oo or x — > + ; for 
this to occur, the timescales of such dynamical instabilities 
should be shorter than the timescale of self-similar evolu- 
tion around Ti. Such dynamical instabilities can give rise to 
a neutron star. This might then be viewed as a two-stage 
process to produce a neutron star with a progenitor mass 
roughly in the range of ~ 6 — 8Mq or higher. 

For progenitor stars with masses less than ~ 6 — 
8Mq, similar rebound shocks may occur but with weaker 



strengths. Such rebound shocks could fail to do the driv- 
ing completely as a result of various energy losses, yet their 
presence and signatures may be detectable. 

We now explore a few more examples for producing 
proto-white dwarfs or proto-neutron stars in core collapse 
and rebound shock processes. For all these cases below, we 
now take rj — 3 x 10 8 cm and r = 10 13 cm. 

The relevant parameters for the first rebound shock 
model are: n = 0.7 (or equivalently 7 = 1.3), ao = 0.175, 
xq = 1.672, wo = 0.293, x S2 = 1.1, A = 1.133, A = 0.752, 
B = -0.222, k 2 = 1.65 x 10 17 cgs unit, ki = 1.29 x 10 17 
cgs unit, M 0i i n i = 7.85-Mq, Mi jU u — 1.36M0. This is a case 
of producing a remnant proto-white dwarf near the Chan- 
drasekhar mass limit of Mi tU lt = 1.36Mq from a progenitor 
star of mass M 0t i„i = 7.85Mq. We feel that with more phys- 
ical input and requirements, our model can provide sensible 
dynamic constraints on the initial and final masses of pro- 
genitor and of remnant proto-white dwarf. 

The relevant parameters for the second rebound shock 
model are: k 2 = 2.2 x 10 17 cgs unit, k\ = 1.71 x 10 17 cgs 
unit, M ,ini = H.8M0, Mi t ult — 2.O5M0; other parameters 
are the same as the first model parameters above. As al- 
ready noted earlier, this proto-white dwarf may be subject 
to dynamical instabilities and could lead to the formation of 
a proto-neutron star in the end. 

The relevant parameters for the third model are: n = 
0.7 (or equivalently 7 = 1.3), ao = 0.196, xo — 1.324, Vo = 
0.0341, x 32 = 0.8, A = 1.146, A = 0.5585, B = -0.864, 
k 2 = 2.1 x 10 17 cgs unit, ki = 1.60 x 10 17 cgs unit, M ,im = 
7.96M Q , M itU it = 1.92M Q . Again, this proto-white dwarf 
may be subject to dynamical instabilities and could lead to 
the formation of a proto-neutron star eventually. 



9 DISCUSSION AND SPECULATIONS 

We are able to construct various semi-complete self-similar 
solutions either with or without a rebound shock; these solu- 
tions may be adaptive to various astrophysical settings. Our 
numerical examples may be interpreted as an approximate 
model of a rebound shock initiated during the gravitational 
core collapse in a massive star, which drive materials out to 
eventually form possible quasi-static configurations. 

For the convenience and clarity of discussion, we mainly 
focus on rebound shocks for supernovae and remnant neu- 
tron stars in the preceding sections. In fact, we may well 
choose a larger value of k 2 for shock models 3 and 4 in Figure 
|H| say, 5 times larger, it follows that all the masses and ra- 
dial flow speeds should be multiplied by the same factor of 5. 
Then the initial progenitor mass for these two shock models 
would be around 35M0 and the final mass enclosed within 
n would be 8.25M0. This then implies a black hole left 
behind as the remnant compact object during the core col- 
lapse and rebound shock processes 6 . In addition, we can also 
choose relevant parameters to model the formation of white 
dwarfs, for instance, taking n — 3 x 10 8 cm and r = 10 13 
cm, we can construct an example for M 0i i„i — 6.3M0 and 
Mi tU u = I.IM0. Although we provide quantitative exam- 
ples using our rebound shock model framework, our main 

e Brown & Bethe 1994 suggested supernova produced low-mass 
black holes with masses only slightly above 1.5Mq. 
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motivation is to point out various conceptual possibilities 
seemingly physically plausible. After all, our current hydro- 
dynamic model is very limited in many ways. 

While we cannot draw sure conclusions on the fate 
of progenitor stars, we may outline a scheme for possible 
remnant compact objects left behind based on the current 
theoretical knowledge. Within the inner reference radius 
Ti ~ 10 6 cm, the final static configuration is assumed to be 
a proto-neutron star if the enclosed mass is in the range of 
< 3 — 4Mq or a black hole if the enclosed mass is higher 
than ~ by the same token, within the inner reference 

radius n ~ 3 — 4 x 10 8 cm, the final static configuration is 
assumed to be a proto-white dwarf if the enclosed mass is 
in the range of < 1.4M© or a proto-neutron star if the 
enclosed mass is higher ~ 1.4Mq. 

Our analysis on the sonic critical curve shows that the 
critical curve has qualitative differences as compared to that 
of the isothermal case (e.g., Shen & Lou 2004), i.e., at small x 
values the sonic critical curve diverges here, and there exists 
a lower limit of the reduced density a on the sonic critical 
curve. Such analysis is technically useful in determining the 
sonic critical curve, and is also adaptive in determining the 
asymptotic behaviours of other complicated functions. 

The 'quasi-static' asymptotic solution is so named af- 
ter the first-order of this solution (viz., the semi-complete 
global SPS solution 1111 or l2Tfll . This new asymptotic solu- 
tion at small a; is a characteristic of polytropic flows only. 
For the first type of this 'quasi-static' asymptotic solution 
with a real K > 1, since L can be either positive or neg- 
ative, the perturbation on the SPS can be either outflow 
or inflow and influences the density profile correspondingly 
(density perturbation is positive for an outflow and nega- 
tive for an inflow). The second type of this 'quasi-static' 
asymptotic solution with complex K represents a wave-like 
perturbation imposed on a SPS. Both the radial flow speed 
and density profiles have an oscillatory feature, and accord- 
ing to equations 1341 and 1351 1 . the oscillatory terms have 
a phase difference not larger than 7r/4 between density and 
velocity profiles. The x Kl exp(iK2 In x) term in radial flow 
speed profile has an interesting property that the vibration 
decreases when x gets smaller. 

In the context of supernova explosions, we propose to 
search for characteristic signatures in density profiles of both 
types of such 'quasi-static' asymptotic solutions in the Cas- 
siopeia A supernova remnant as well as other suitable su- 
pernova remnants. While unconventional, we suspect that 
some proto-white dwarfs are also formed involving rebound 
shocks and explosions. 

Since various other physical conditions, such as the neu- 
trino opacity, radiation pressure, general relativistic effects, 
magnetic field, and rotational effects are generally involved 
in such core collapsing phase of stellar evolution, we hope 
that this scenario and interpretation, while limited by a 
highly idealized model framework, catch some essential fea- 
tures of supernova explosion. This scenario should be tested 
by observations and numerical simulations. We plan to in- 
corporate a random magnetic field into this model frame- 
work. Thus a pa rallel analysis of a magnetohydrodynamic 
model similar to lYu fc Loul (I2005T) and Yu et al. (2006) is 
anticipated. 
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„ 2a! 2(a') 2 1 2/s 3 k/o _ , . _. 

a + h—^ - -a 2/i + -a o/i = . (A5) 

x 3a 2 4 

We shall pursue a more complete analysis of this problem in 
a separate paper. 



APPENDIX A: POLYTROPIC CASES OF 7 > 4/3 

When index 7 > 4/3 is adopted for a polytropic gas, the 
self-similar transformation should be modified slightly. In 
fact, the only difference between this case and what we have 
considered in the main text concerns the transformation for 
the total enclosed mass M. As M > is a physical require- 
ment in transformation (7J, the scaling factor e should take 
the form of 

6 = (2^3njG (A1) 

for 7 > 4/3 in a parallel analysis. For example, the nx~v > 
criterion for m > should now be replaced by nx — v < 0. 

The 7 = 4/3 case is special and makes the scaling factor 
e in definition J7J meaningless. Noting that the factor 3n — 2 
was only added for convenience, we may simply set 

1.3/2 

to perform a self-similar transformation. Substituting this 
modified transformation into hydrodynamic equations Q 
through under spherical symmetry, we obtain 

v = nx = 2x/3 , (A3) 

m = ax 2 , (A4) 



and a second-order nonlinear ODE for a alone 



